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Consider a set of A'^ systems and an arbitrary interaction Hamiltonian H that couples tiiem. 
We investigate the use of local operations and classical communication (LOCC), together with 
the Hamiltonian H, to simulate a unitary evolution of the A'' systems according to some other 
Hamiltonian H' . First, we show that the most general simulation using H and LOCC can be also 
achieved, with the same time efficiency, by just interspersing the evolution of H with local unitary 
manipulations of each system and a corresponding local ancilla (in a so-called LU+anc protocol). 
Thus, the ability to make local measurements and to communicate classical information does not 
help in non-local Hamiltonian simulation. Second, we show that both for the case of two d-level 
systems (d > 2), or for that of a setting with more than two systems (A*' > 2), LU+anc protocols 
are more powerful than LU protocols. Therefore local ancillas are a useful resource for non-local 
Hamiltonian simulation. Third, we use results of majorization theory to explicitly solve the problem 
of optimal simulation of two-qubit Hamiltonians using LU (equivalently, LU+anc, LO or LOCC). 

PACS numbers: 03.67.-a, 03.65.Bz, 03.65.Ca, 03.67.Hk 



I. INTRODUCTION 

The problem of using a given non-local Hamiltonian 
H and some class of local operations to simulate another 
non-local Hamiltonian H' has very recently attracted the 
attention of several authors in quantum information sci- 
ence ^, H, ^, 1^ . Nonetheless, average Hamiltonian 
techniques, a basic ingredient in non-local Hamiltonian 
simulation, have been studied for many years in control 
theory , and are commonly used in the area of nuclear 
magnetic resonance 

From the perspective of quantum information science, 
non-local Hamiltonian simulation sets a frame for the 
parameterization of the non-local resources contained in 
multi-particle Hamiltonians, very much in the line of 
thought pursued to quantify the entanglement of multi- 
particle quantum states. In the most common setting, 
fast local unitary operations LU are performed on a se- 
ries of systems to effectively modify the Hamiltonian H 
that couples them. A remarkable result is the qualita- 
tive equivalence of all bipartite interactions under LU 
[||, ^ |, H, 0]. This can be shown to imply that any 
Hamiltonian H with pairwise interactions between some 
of the systems can simulate any other Hamiltonian H' 
consisting of arbitrary pairwise interactions between the 
same systems. 

At a quantitative level, the time-efficiency with which 
a Hamiltonian H is able to simulate a Hamiltonian H' 
can be used as a criterion to endow the set of non-local 
Hamiltonians with a (pseudo) partial order structure, 
that allows to compare the non-local capabilities of H 
and H' jj]. For two-qubit Hamiltonians, simulations us- 
ing LU or arbitrary local operations LO have been shown 
to yield the same optimal time efficiencies, and the result- 
ing partial order structure has been computed explicitly. 
This has led to the necessary and sufficient conditions for 



H to be able to simulate H' efficiently for infinitesimal 
times, that is, the conditions under which the use of H 
for time t allows to simulate H' for the same time t, in 
the small time t limit. Equivalently, this result shows 
how to time- optimally simulate H' with H, in the sense 
of achieving the maximal simulation ratio t' /t, where t is 
the time of interaction H that it takes to simulate inter- 
action H' for a time t' . 



A. Ancillary systems, generalized local 
measurements and classical communication in 
non-local Hamiltonian simulation 

The aim of this paper is to elucidate the role a number 
of resources play in the simulation of non-local Hamil- 
tonians. Relatedly, we seek at establishing equivalences 
between different classes of operations that may be used 
in a simulation protocol. 

We first address the question whether classical commu- 
nication (CC) between the systems is useful in non-local 
Hamiltonian simulation. Recall that in protocols that 
include local measurements, the ability to communicate 
which outcome has been obtained in measuring one of 
the systems allows for subsequent operations on other 
systems to depend on this information. Now, can this 
ability be used in non-local Hamiltonian simulation to 
enlarge the set of achievable simulations? Suggestively 
enough, the answer is yes in the closely related problem 
of converting one non-local gate into another non-local 
gate using LO. For instance, a series of two-qubit gates 
U exist such that they can be achieved by performing a 
C-NOT gate and LOCC but can not be achieved by a 
C-NOT and LO @. 

We also study the advantage of using ancillary systems 
in simulation protocols, as well as performing general 
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local operations instead of just local unitary transfor- 
mations. Altogether, our analysis refer to the following 
classes of transformations: 

• LU = local unitary operations, 

• LU+anc — local unitary operations with ancillas, 

• LO = local operations 

• LOCC ~ local operations with classical communi- 
cation. 



B. Results 

This paper contains the following three main results 
concerning the simulation of non-local Hamiltonian evo- 
lutions for infinitesimal times: 

(i) LOCC (or LO) simulation protocols can be reduced 
to LU-f anc simulation protocols. That is, for A^-particle 
Hamiltonian interactions H and iJ', any protocol that 
simulates H' using H and LOCC (or LO) can be replaced, 
without changing its time efficiency, with a protocol in- 
volving only H and local unitary transformations. Each 
local unitary transformation may be performed jointly 
on one of the TV systems and a local ancilla. 

(m) Apart from exceptional cases such as that of two- 
qubit Hamiltonians Q — in which any LU-|-anc protocol 
can be further replaced with an even simpler protocol 
that uses only LU on each qubit — , the use of ancillas is, 
in general, advantageous. This is proven by constructing 
explicit examples of LU-f-anc protocols where ancillas are 
used to obtain simulations that cannot be achieved with 
only LU operations, both in the case of two d- level sys- 
tems (d > 2) and in the case of > 2 systems. 

(iii) For two-qubit Hamiltonians, we use results of ma- 
jorization theory to recover the optimality results pre- 
sented in In view of the equivalence between LU, 
LU-|-anc, LO and LOCC protocols for two-qubit sys- 
tems, this solves the problem of time-optimal, two-qubit 
Hamiltonian simulation under any of these classes of op- 
erations. 

The structure of the paper is as follows. In section II 
we introduce some known results. Section III, IV and 

V present results (i), (ii) and {iii) respectively. Section 

VI contains some conclusions and appendices A and B 
discuss some technical aspects of sections III and V. 



II. PRELIMINARIES 

We start by reviewing some background material from 
Ref. 1^, of which the present work can be regarded as 
an extension. 



A. Non-local Hamiltonian simulation and classes 
of operations 

Recall that the aim of non-local Hamiltonian simula- 
tion is, given a set of systems that interact according to 
Hamiltonian H for time t and a class C of local control 
operations, to be able to produce an evolution e~*^ * for 
the systems, where H' and t' are the simulated Hamilto- 
nian and the simulated time. [We take h = 1 along the 
paper] . 

As mentioned above, one can consider several classes 
of operations to assist in the simulation, including LU, 
LU-|-anc, LO and LOCC. As in [Q, we make two ba- 
sic assumptions: (i) these additional operations can be 
implemented very fast compared to the time scale of the 
Hamiltonian H (we actually consider the setting in which 
they can be performed instantaneously and thus charac- 
terize the fast control limit); (ii) these operations are a 
cheap resource, so that optimality over simulation pro- 
tocols is defined only in terms of the ratio t' /t, that is, 
in terms of how much time t' of evolution according to 
H' can be produced by using H for a time t. Another 
interesting parameter characterizing simulations, that we 
do not analyze here, would be some measure of the com- 
plexity of the simulation, that is of the number of control 
operations that are performed. 

We also note that the inclusions between classes of op- 
erations, LU C LU-|-anc C LO C LOCC, imply relations 
between the sets of achievable simulations and time effi- 
ciencies. For instance, since LOCC simulation protocols 
contain all LU simulation protocols, we expect LOCC 
protocols to be more powerful than LU protocols. 

B. Infinitesimal— time simulations 

The maximal simulation factor s{t') = t' /t when sim- 
ulating e~*^ * by using H for time t may depend on t' . 
However, we are ultimately interested in characterizing 
the non-local properties of interaction Hamiltonians, irre- 
spective of interaction times. A sensible way to proceed 
is by considering the worst case situation, namely the 
time t' for which the optimal ratio s{t') achieves its min- 
imal value. This occurs for an infinitesimal time t' . That 
is, simulations of H' for a time such that << 1 

are, comparatively, the most expensive in terms of the 
required time t of interaction H . The reason is that, (i) 
simulations for an infinitesimal time are a particular case 
of simulation, providing an upper bound for the mini- 
mum of s(t'), and {ii) any finite-time simulation — or 
gate synthesis — can be achieved, maybe not optimally, 
by concatenating infinitesimal-time simulations. 

We shall denote Sh'\h the limit limt'^o s{t'), and call 
it the simulation factor of H' with H [sh'\h corresponds 
to the inverse of the time overhead fi of Ref. Q, that 
is Sff/^fj = Then, apart from quantifying the time 

efficiency in infinitesimal simulations, s^f'i^f has also two 
other meanings: 



3 



T' /sh'\h upper bounds the time T of use of H 
needed to perform the unitary gate e~*_^^ , for any 



T' {gate simulation or gate synthesis ,15 1); 



sh'\h is the optimal time efficiency in dynamics 
simulation. That is, sh'|h is the maximal achiev- 
able ratio T' /T , where T is the time of H required 
to simulate the entire evolution of a system accord- 



ing to e 



-it'H' 



where t' runs from to T' |l^ 



In an abuse of notation, we shall refer to condition 
ll-ff't'll << 1 as the smah time limit, of which 0{t') [or 
0{t)] will denote first order corrections. 



C. Optimal and efficient simulations 

For any class C e {LU, LU+anc, LO, LOCC} of the 
above operations and in the small time limit, the space 
of achievable evolutions using Hamiltonian H and oper- 
ations C turns out to be convex. Then the following two 
problems, 

PI: Given any H and H' , determine when H' can be 
efficiently (i.e., t' = t) simulated with H for infinitesimal 
times, denoted 



III. EQUIVALENCE OF LOCC AND LU+ANC 
PROTOCOLS 



In this section we show that any protocol for non-local 
Hamiltonian simulation based on LOCC can be replaced 
with another one based only on LU-|-anc and having the 
same time efficiency. This result, valid for infinitesimal- 
time simulations on arbitrary A^-particle systems, brings 
an important simplification to the general problem of 
non-local Hamiltonian simulation, since it implies the 
equivalence of LOCC, LO and LU-fanc protocols. 

We first describe in detail a most general protocol for 
Hamiltonian simulation using LOCC. Then we show — 
through an argument that exploits the fact that entangle- 
ment only decreases under LOCC — that any such proto- 
col can be replaced with another one using only LU-|-anc. 
The key point of the proof is to assume that one of the 
systems is initially entangled with an auxiliary system 
Z, and to realize that a non-trivial measurement (i.e., a 
measurement not equivalent to some local unitary trans- 
formation) on the system would partially destroy this 
entanglement in an irreversible way. Since we are sim- 
ulating a unitary process on the systems (which should 
preserve the entanglement between those and Z), all lo- 
cal measurements must be trivial, and can be replaced 
with unitary transformations. 



H' >r H: 



(1) 



P2: Given any H and H' , determine the simulation fac- 
tor Sh'\h; 

are equivalent, since sh'\h is nothing but the greatest s 
such that sH' can be efficiently simulated by H, that is, 
such that sH' >c H. 

D. Equivalence of LO and LU+anc protocols 



A. Hamiltonian simulation using LOCC 

For clearness sake we will perform most of the analysis 
in the simplest non-trivial case, that involving only two 
qubits, because this already contains all the ingredients 
of the general A^-particle setting. Let us consider, then, 
that qubits A and i3, with Hilbert spaces TLa and TYb, 
interact according to H for an overall time t, and that, 
simultaneously, they are being manipulated locally. 



The simulation of non-local Hamiltonians using LO 
and that using LU+anc are equivalent (see [Q for de- 
tails) , in that any protocol based on LO can be replaced 
with another one that uses only LU+anc and that has 
the same time efficiency. The ultimate reason for this 
equivalence is that even if LO provide, through measure- 
ment outcomes, information that can be used to decide 
on posterior local manipulations, this information cannot 
be transmitted to the other parties [unless the interaction 
itself is used for this purpose, but this leads to null effi- 
ciency t' /t when t — > 0]; then, unitarity of the simulated 
evolution implies that each party is effectively applying 
a trace-preserving local operation on its subsystem, and 
this can always be achieved using only LU+anc. 

The previous situation changes when classical commu- 
nication is allowed between the parties, because then they 
can coordinate their manipulations. In spite of this fact, 
CC does not help in Hamiltonian simulation, as we move 
to discuss next. 



1. Local manipulation 

The most general local operation on, say, qubit A can 
be achieved by (i) appending to A an ancillary sys- 
tem A' in some blank state [0^') S Ti-A'] (**) performing 
a unitary transformation U on TLaa' = "Ha ® T^A'] {Hi) 
performing an orthogonal measurement on a factor space 

rneas of the total Hllbcrt space T~CaA' — /C T'^meas^ 

given by projection operators {P°'}] and {iv) tracing out 
a factor space T" of Haa' — 'Hout ® where Tig^i 
and T" may depend on the measurement outcome a of 
step {iii). Under {i)-{iv) the initial state \4>a) of qubit 
A transforms with probability pa according to 

|0a)(<^a| --fa(|0A)(<^A|) = 

— trr^ [P"C/(|0a)(0a| <^ |0a.)(0a'|) C/^P"] , (2) 

Pa 
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where p„= tr[P"f7(|,^^)(0^| ® |0a')(0a' We 
can introduce operators M" : Ha T^out^ 

Aff = (i"|F"L/|Oa'), (3) 

where {!«")} is an orthonormal basis of T". Then Eq. 
can be rewritten as 

U\^a){M) ^-J^MnMiMM^^. (4) 

Pa . 

Now, since in our case the eventual result of this manip- 
ulation must be a unitary evolution, we are interested in 
transformations £a that map pure states into pure states, 
that is, such that can be implemented by just one oper- 
ator M° : Ha -> H'^uu 

\<Pa){(^a\ ^ —A'ncj>A){(I^A\M^\ (5) 

Pa 

Pa = tr[M'^\(j)A){(l)A\M"^]. Therefore the effect of the 
local manipulation on qubit A is a generalized measure- 
ment M that, with probability p", maps the state of A 
into a state supported on 7i";^j, 

\cj>A) -^M'^\^a), (6) 

and produces classical information a. The measurement 
operators {Af"} characterizing M satisfy Af^tM" = 

I A. 

More generally, in a simulation protocol measurement 
M may depend on some previous information (3, in which 
case we write M.^ . In addition, the corresponding mea- 
surement operators {Af "'^} may map states from a two- 
dimensional subspace 7i,f„ C Haa' into another two di- 
mensional subspace 7^"^^ C Haa' that depends both on 
the measurement outcome a and on the previous infor- 
mation /3, that is 

M'^^^ : Ff„ ^ Ki- (7) 

In the following, a series of measurements M. will be 
concatenated, in such a way that the out-subspace Hout 
for a given measurement is related to the in-subspace Hm 
for the next one. 

We consider that a sufficiently large ancillary system 
A' in a pure state has been initially appended to qubit A 
so that it provides at once the extra degrees of freedom 
needed to perform all generalized measurements M on A. 
Finally, all the above considerations apply also to qubit 
B, to which an ancillary system B' is appended. 

S. LOCC simulation protocol 

A LOCC protocol for simulating e~'* ^ hy H for time 
t is characterize by a partition {ti,t2, ■■■,tn} of t, where 
t'i > 0, ti = t, and a series of local measurements, 
{{Mo,Afo), (Xi\-^r')' iM^-,Af^")}. The protocol 
runs as follows: 



1. The simulation begins with measurements A^o and 
Ao being performed on A and B, respectively. 
These map the original state of AB into a state 
supported on some subspace of AA'BB' . 

2. Then the two qubits A and B are left evolve ac- 
cording to H for a time ti. 

3. After that, measurements A^"^ and Af^^ are per- 
formed. Here, index ai indicates that the measure- 
ments being performed after time ti may depend on 
the outcomes of measurements Mq and Aq. 

4. Again, the measurements are followed by an evolu- 
tion, for time t2, of yl and B according to H, and 
the protocol continues in an iterative fashion. 

5. In step k, qubits A and B are first left evolve ac- 
cording to H for a time tk and then measurements 
A^J!* and Af^'' (a^ denoting again a possible de- 
pendence on the outcome of any previous measure- 
ment) are locally performed in AA' and BB'. 

6. The protocol finishes after measurements A^^" and 
A/""" have been performed. These last measure- 
ments must leave the two-qubit system AB in a 
pure state (that is, uncorrelated from systems A'B', 
that are traced out). 

Thus, the two-qubit system AB is initially in some 
state becomes entangled with the ancillas A' and B' 
during the manipulations describes above, but ends up 
in the state e^'^ * after time t. 

Note that the protocol described above has a tree 
structure, starting with a preestablished couple of local 
manipulations and ending up at the extreme of a branch 
characterized by the outcomes of all (conditioned) lo- 
cal operations performed during the time interval t. We 
move now to characterize one of these branches. 



3. One branch of the protocol 

Let us suppose we run the simulation once. This cor- 
responds to some given branch of the protocol, that we 
label r, and have represented in the figure. Branch 
F is characterized by a series of measurement opera- 
tors {{M^,N^),...,{M^,N^)}, where the superindices 
ckfc containing the information that characterizes the 
branch have been replaced with F to simplify the nota- 
tion. Recall that the aim of the protocol is to achieve an 
evolution according to e~'^ * . Therefore, for any initial 
vector 1^) of the two-qubit system AB, the measurement 
operators {{M^ , Nl)}^^^ must obey 

VpFe-'^'*' \^) = {Ml ® Oe-""^ • • • 

• • • [Ml ® Nl)e-''^"{M^ ® Nl)\^j), (8) 

where py denotes the probability that branch F occurs in 
the protocol. Eq. (^ is the starting point for the rest of 
the analysis in this section. 




FIG. 1: Schematic representation of a Hamiltonian simulation protocol using LOCC. The unitary evolution of the composite 
system AB according to H and for a time t = ti is interspersed with local measurements TM^*" [on systems AA'] and Af^'' [on 
systems BB'] to obtain a unitary evolution of AB according to H' and for a time t' . Here at indicates the local measurements 
performed at step k, which may depend on the outcome of all previous steps thanks to the ability to communicate classical 
between the systems (winding lines). In the figure we have replaced the superscripts Ofc with F, F denoting a particular branch 
of the protocol (cf. Eq. (^)). Thus, in branch F measurement operators [corresponding to measurement M]^] and A'^^ 
[corresponding to measurement Af^] transform systems AA' and BB' at step k of the protocol. 



B. LOCC protocols are as efficient as LU+anc 
protocols for infinitesimal time simulations 

As discussed in the introduction, we are interested here 
in simulations for an infinitesimal simulation time t. In 
this regime Eq. (|^) significantly simplifies, because we 
can expand the exponentials to first order in t (or equiv- 
alently, in {tk} and t'), thereby obtaining an equation 
which is linear both in H and H' . In addition, if t is 
small then qubits A and B interact only "a little bit". 
In what follows we will use this fact to prove the main 
result of this section, namely that all the measurement 
operators {M^, N^}^^q in Eq. (||) must be, up to neg- 
ligible corrections, proportional to unitary operators in 
some corresponding relevant supports. This will eventu- 
ally imply that LU+anc protocols can already simulate 
any evolution e~*^ * achievable in a LOCC protocol. 

We note that this result is not valid for the intercon- 
version of non-local gates jl^. There the systems are 
allowed to interact according to a finite gate (e.g. a C- 
NOT), and thus accumulate some finite amount of en- 
tanglement (e.g. an ebit) in the ancillary systems, that 
can be used, together with LOCC, to perform some new 
non-local gate (e.g. through some teleportation scheme). 



1. LOCC protocols for infinitesimal-time simulations 
We define a series of operators Mk and Af^ by 



Ml. = Mi • • • M, 



r 

k ' 



k = l. 



Mo = Ml---Ml, 



and also an analogous series of operators Nk,N'^ and iVo. 
Notice that operator M'f. describes a concatenation of all 
local measurements in branch T performed from the be- 
ginning of the protocol and up to step fc — 1 on the state 
initially supported on Ti-AA'^ while Mk collects the ma- 
nipulations that will be performed from step k until the 
end of the protocol. In the small time regime, we can 
expand the exponentials in Eq. (^) as a series in tk and 
t' to obtain, up to second order corrections C'(i^), 



^{lAB-istH')\^) ^ (10) 
Mo ® iVo - itJ2Pk{Mk <E) Nk)H{M'k ® K)] 1^), 



k=l 



where we have introduced probabilities pk = tk /t and the 
efficiency factor s = t' /t oi the branch, so that all times 
are expressed in terms of t. 
This equation indicates that 



Mo ® iVol^) = ^lABlij) + 0{t), 



(11) 



for any two-qubit state ji/)), from which it follows that 
the probability pr that branch T occurs cannot depend 
on 1-0) up to 0{t) corrections, also that both Mq and A^o 
must be proportional to the identity operator in Ha and 
Hb, 



Mq - ^qlA + Oit), 
No = q-'lB+0{t), 



(12) 
(13) 



(9) 



where q is some positive parameter. Notice that the order 
t corrections in Eq. (11) correspond to local terms, that 
is, to operators of the form t{lA ® Ob + <8) Ib), and 
thus are irrelevant to this discussion ||T5|| . In what follows 
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we neglect these local terms for clearness sake. Bea ring 
this remark and Eq. (^l|) in mind, we rewrite Eq. ( |lO| ) 
as the operator equation 



{Iab - istH') ^ (14) 
P^Iab ~ itY,Pk{Mk ® Nk)H{Ml ® iV^) + 0(1^). 

k=l 



That is, 



P^sH' = Y^PkiMk ® Nt,)H{M'k K) + 0{t), (15) 
fc=i 



where, because of Eq. (|1 1|) , some other constrains apply. 
More precisely, if and N^. are given by 



Ml = ^q{\^xl){QA\ + \^i\){lA\) 
q-\\4){0B\ + W'l){lB\). 



(16) 



where {I^a)} and are orthonormal basis of Ha 

and Hb and {\iJl) € Haa'} and {\vi) G TIbb'} are ar- 
bitrary vectors, not necessarily normalized, then Mk and 
A^fc must fulfill 



A'h = |0a)(mSI + |1a)(Ai1+OW, 
Nk = \OB){i^'o\ + \U){i>',\ + 0{t), 



(17) 



where {I/if )} is the biorthonormal basis of {I/if )} (in the 
subspace spanned by {|m?)})7 that is (/i*^|/t*^) = Sij, and 
similarly {pf}} is the biorthonormal basis of {|i^}, so 
that Mo^No = {MkM'^) ® [NkN'^) fulfills Eq. (y). 

Now, going back to the measurement operators Af^, 
we can expand them as 

Ml = ^q{\^Jil){OA\ + \^A){U\ ) 

Ml = |/^r')(/iSl + lA^^')(Mj| A; = l,...,n-1 
Ml = \QA){m + \^A){m+0{t) (18) 

and similarly for the Nl . 

2. Unitarity and conservation of entanglement 

We carry on this analysis by focusing our attention 
only on the operations performed on systems AA' . We 
will show that operators Mk and Mj, can be replaced with 
operators proportional to (O^'l^fe and uI\Oa'), where Uk 
is a unitary matrix acting on Haa' ■ We will use the fact 
that the protocol must be able to keep the entanglement 
of A with another system Z. 

Let us suppose, then, that qubit A is entangled with a 
distant qubit Z , with the maximally entangled vector 



V2 



(|0a)®|0z) + |1a)® |lz)) 



(19) 



describing the pure state of AZ. Any unitary evolution 
of qubits A and B preserves the amount of entanglement 



between qubit Z and qubits AB. In particular, if the 
unitary evolutions according to H are infinitesimal, then 
up to 0{t) corrections qubit Z must be still in a max- 
imally entangled state with A after the simulated cvo- 



lution 



This sets very strong restrictions on the 



kind of measurements that can be performed on A during 
the simulation protocol. If during the fc*'* measurement 
in branch T part of the entanglement is destroyed, then 
the simulation protocol necessarily fails with some prob- 
ability, because the destroyed entanglement can not be 
deterministically recovered. Indeed, even if subsequent 
measurement operators in branch T would be able to re- 
store the entanglement and so obey Eq. (^), another 
branch T' diverging from T after the fc*'* must necessar- 
ily fail to recover the entanglement (recall the monotoni- 
cally decreasing character of entanglement under LOCC, 
see e.g. jl^) and thus with some probability the protocol 
must fail to simulate the unitary evolution [ p7[ . 

Let us see the effect of this restrictions on the first 
measurement operator in Eq. (^8|). It transforms 
the initial entangled state into a new state proportional 
to 



ImI) 



II; 



(20) 



which remains maximally entangled if and only if 
|||/io)|| = |||/ii)|| = ri and (/io|/^i) — 0. But this are pre- 
cisely the conditions for M[{— M^) to be proportional 
to a unitary operator from Ti.A to the out-space Tig^j 
spanned by {|/ij^)} or, equivalently, to an isometry from 
Ha to Haa' ■ Thus, we can write 



M[^nUl\OA'), 



(21) 



where 



uI^\^{OaOa' 
r\ 

1=1 



ImI) 



(IaOa' 



ri 

\^i,i){^aIa' 



is some unitary operation defined on Haa'- Here c^a' 

is the dimension of Ha' and {|C;,o), lCi,i)}f=i is a set of 
irrelevant vectors that together with l/tj)/?"! and |/tj)/ri 
form an orthonormal basis of Haa'- Eq. ( [l^ ) implies 
that, in addition, 



Ml = 



prq 



{Oa'\Ui 



(22) 



This characterization in terms of a unitary transfor- 
mation can now be easily extended to the rest of opera- 
tors Mk and Mj,. We use induction over k. We already 
have that the characterization works for fc = 1 . Suppose 
it works for some k — I, that is, in the decomposition 
Eq. ([16|) for M^_i we have |||/tS"^)ll = IIIm?"^) 
{^J'Q^^\^J■l^^) = 0. This means that after the (fc - 
measurement in branch F, the initial state of Eq, 
becomes a state proportional to 



and 

th 



|0. 



o(t). 



(23) 
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where the 0{t) corrections are due to evolutions of AB 
according to H for a time of order t, which slightly en- 
tangle B with AZ . Then, preservation of entanglement 
during the A:*'' measurement (implemented by operator 
Affc_i) requires that also = = and 

(A*fclMfe) = 0; E^^^d therefore 

= nUl\0A')+O{t), 

Mk = y^{0A'\Uk + Oit), (24) 

for some unitary transformation Uk acting on Ti.AA' ■ 

The same argument leads to expressing the operators 
Nk and in terms of unitary transformations Vk acting 
on H-BB' as 

K = skV^\OB') + o{t), 

Nk = —{0B'\Vk+Oit). (25) 
rkq 

Therefore Eq. ( p^ ) finally reads, up to 0{t) corrections 
that vanish in the t — + or fast control limit, 

sH' ^Y.P'^(^A'OBMUk®Vk)H{ul®V^)\OA'OB'). (26) 

k 

3. Equivalence between LOCC and LU+anc protocols 

The set Sjf"~^^ of non-local Hamiltonians that can be 
efficiently simulated by H and LOCC is convex: if H can 
efficiently simulate Hi and H2, then it can also efficiently 
simulate the Hamiltonian pHi + (1 — p)H2. Indeed, we 
just need to divide the infinitesimal time t into two parts 
and simulate Hi for time pt and then H2 for time {l—p)t. 
The resulting Hamiltonian is precisely the above average 
of Hi and H2. Thus, in order to characterize the convex 
set S^occ^ 

we can focus on its extreme points. Notice 
that the previous convexity argument also holds for the 
set S^^"''^'^ of Hamiltonians that can be efficiently sim- 
ulated with LU-|-anc, so that S^+"-'^<^ is also convex. 
Recall also that 3^^+"'"' c S^'^^^ . 

Now, Eq. (|2^) says that all points in Sjf"^''^' can be 
obtained as a convex combination of terms of the form 

{Oa'Ob'HU ® V)HiU^ ® V^)\0a'0b'). (27) 

In addition, in appendix A we show that any such a term 
can be obtained in a simulation protocol using LU-|-anc. 
It follows that (i) any extreme point of Sjf^'~^ is of the 
form (pTf), and that (ii) any extreme point of Sj^'^^'~'' 
belongs to 5'^'^+''"^ so that 5^^+""^ = Sjf^^. This 
finishes the proof of the fact that infinitesimal time sim- 
ulations using LOCC can always be accomplished using 
LU-|-anc. 

Summarizing, we have seen that any (rescaled) two- 
qubit Hamiltonian sH' achievable in branch F of our 
LOCC-simulation protocol (cf. Eq. ([26[)) can also be 



achieved, with the same time efficiency, by just using 
local unitary transformations and ancillas as extra re- 
sources. It is now straightforward to generalize the above 
argument to N systems, each one having two or more 
levels, thereby extending the equivalence of LOCC and 
LU-|-anc protocols to general multiparticle interactions. 
Indeed, for any d-level system involved in the simula- 
tion, we just need to require that its entanglement with 
some remote, auxiliary d-level system be preserved, and 
we readily obtain that all measurements performed dur- 
ing the simulation protocol can be replaced with local 
unitary operations. We thus can conclude, using the no- 
tation introduced in section II. B, that 

H' >LOCC H <^ H' > LU+anc H. (28) 



IV. LU+ANC PROTOCOLS ARE NOT 
EQUIVALENT TO LU PROTOCOLS 

The equivalence between infinitesimal-time simula- 
tions using LOCC and LU-f-anc may be conceived as 
a satisfactory result. On the one hand, it discards lo- 
cal measurements and classical communication as useful 
resources for the simulation of non-local Hamiltonians. 
This essentially says that in order to simulate Hamilto- 
nian dynamics, we can restrict the external manipula- 
tion to unitary operations, possibly involving some an- 
cillary system. In this way the set of interesting sim- 
ulation protocols has been significantly simplified. On 
the other hand, it is reassuring to see that, despite the 
diversity of classes of operations that we may use as 
a criterion to characterize the non-local properties of 
multiparticle interactions, most of these criteria (LOCC, 
LO and LU-|-anc) yield an equivalent classification and 
quantification. In other words, we do not have to deal 
with a large number of alternative characterizations. We 
shall show here, however, that simulation using only LU, 
that is, without ancillas, is not equivalent to that using 
LU-l-anc. 

The reason for this inequivalence is the following. Con- 
sider a multi-partite Hamiltonian of the form Ha ® 
Hbc -1 where Ha acts on a d dimensional space T-La and 
Hbc- - acts on TLb®'Hc ■ ■ ■ ■ In the presence of an ancilla 
Ha' , LU can be used so that operator Ha acts on some 
d dimensional factor space IC of Haa' {H,aa' — K. ® K.'). 
The net result is an effective Hamiltonian acting on T-La- 
As the following examples show, some of these effective 
Hamiltonians can not be achieved (at least with the same 
time efficiencies) by using only LU. 



A. LU+anc protocols versus LU protocols 

In the previous section we saw that, in the fast con- 
trol limit, the extreme points of the convex set S]^'^"''^^ 
of bipartite Hamiltonians that can be efficiently simu- 
lated with H using LU+anc, [equivalently, those of the 
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set S];f"-^'-^] are, up to local terms, of the form 

8{H) = {{)a'Qb'\U®V{H ®Ia'b')U'^®V^\'^A'Qb') (29) 

(an analogous expression holds for the multi-partite 
case). Notice that in Eq. ( p9|) we have replaced oper- 
ator H of Eq. ( p7| ) with H (g) Ia' b' to make more explicit 
that ancillas are being used. 

Can all simulations of this type be achieved by us- 
ing only LU? The most general simulation that can be 
achieved from H and by LU reads (see Q for more de- 
tails) 



+ Ib + lA®n + alAB , 



(30) 



where {pk}, X^^-Pfc = 1; '^^ ^ probability distribution, 
{uk\ and {vk\ are local unitaries acting on A and _B, 
m and n are self-adjoint, trace-less operators and a is a 
real constant. The previous question translates then into 
whether for any U and V in Eq. (p9[), we can find a set 
{Pk^'u^k^'Vk\^ 71 and fl such that Eq. ( ^0[ ) equals £{H) 
in (|29|^ 

In^ it was shown that, in the particular case of two- 
qubit systems, the previous conditions can always be ful- 
filled. Next we shall show that this is sometimes not the 
case for Hamiltonians of two d-level systems for c? > 2, 
and also for Hamiltonians of more than two systems. 



B. Inequivalence between LU+anc and LU 
protocols 

1. Example 1: two d-level systems (d> 2) 

We first consider two c?~level systems A and B, d> 2, 
that interact according to 



(31) 



where Pi ® Pj = IiaXt-a] ^ |js)(js|- We will show that 
by means of LU-(-anc, Hamiltonian K can be used to 
efficiently (that is, with unit efficiency factor s) simulate 



K' = PQ^Pl+J2P^®P^■ 



m 



As we discuss in appendix A, the transformation of a 
Hamiltonian H acting on AB, 

£{H) = {OA'OB'PiH ® Ia'B'W^Oa'Ob'), (34) 

can be achieved using LU+anc [notice that this corre- 
sponds to choosing Vbb' = Ibb' in Eq. (Al)]. In par- 
ticular, this transformation takes any term of the form 
Pi (g) Pj into 





E{P^®Pj) = { {Po + Pi)®Pj 

P^ (g Pj 

which in particular implies 

£{K) = K'. 



i = 
i = 1 
i > 1. 



(35) 



(36) 



Now, if this simulation is to be possible with the same 
time efficiency by using only LU, then we must have, 
because of Eq. (|3C|), 



K' = Q + m I® Ib + Ia ® n + gIab , 



(37) 



where Q = Yfi=lHkP^'^kPtu\ g) vuPiv]. > 0, but this 
is not possible. Indeed, we first notice that, by taking 
the trace of this expression we obtain a = 0, whereas by 
tracing out only system B we obtain 

I = I + dm, (38) 

and thus m — 0. Tracing out only system A leads to 



so that n = {—Pq + Pi)/d and condition 

d-l r 



(39) 
becomes 



K' ^ Po®Pi+Y,P^®P^ = Q + -T®{-Po + Pi)- (40) 



Then, recalling positivity of Q, we obtain the following 
contradiction 

= tr [P2 (g PiK'] = tr [P2 ® PiQ] 
+ tr [(P2(gPi) (^^(g(-Po + Pi)^] 
= tr [(Pa (g Pi)Q] + l/d > 1/d. (41) 



We will also show that K' can not be efficiently simulated 
using only LU. 

Let A' be a c?-level ancilla. We need a unitary trans- 
formation U satisfying 



d-l 



{Oa'\U - \Oa){IaOa'\ + 5Z \iA){iAiA'\- (33) 



i=l 



Thus, for any d > 2, we have explicitly constructed 
an example of LU-f anc simulation for Hamiltonians act- 
ing on two d-level systems that can not be achieved using 
only LU. We recall, however, that for two-particle Hamil- 
tonians, LU-|-anc and LU protocols only differ quantita- 
tively, for LU protocols are able to simulate any bipartite 
Hamiltonian H' starting from any other H with non- 
vanishing s^'IH |, § § § 0- 



9 



2. Example 2: a 2 x 2 x 2 composite system 

Let us consider now the simulation, for an infinitesimal 
time t, of the three-qubit Hamiltonian 

'/ 



by the Hamiltonian 



= (73 (g) (73 (g) (73 , 



where 



(73 



1 

-1 



(42) 
(43) 

(44) 



This is possible, when allowing for LU+anc operations, 
by considering the transformation U acting on qubit A 
and on a one-qubit ancilla A' in state lOyi'), where 



{Qa'\U = |0a)(0a| ® (Oa'I + |U)(Oa| ® (U' 

Indeed, we have that (Oa'|C^((73 ® Ia')U^Qa') 
that 

{{)a'\UKU^\Qa')=K'. 



(45) 

I A, SO 

(46) 

On the other hand it is impossible to simulate K' by K 
and LU, for it would imply to transform (73 into / through 
unitary mixing, which is a trace-preserving operation. It 
is straightforward to construct similar examples in higher 
dimensional systems, and also with more than three sys- 
tems. 

We note that, as far as interactions involving more 
than two systems are concerned, the inequivalence be- 
tween LU-|-anc and LU simulation protocols is not only 
quantitative, leading to different simulation factors, but 
also qualitative. The last example above shows that LU 
protocols can not be used to simulate Hamiltonians that 
can be simulated using LU+anc and the same interaction 
H. 



OPTIMAL SIMULATION OF TWO-QUBIT 
HAMILTONIANS USING LOCC 



In this last section we address the problem of optimal 
Hamiltonian simulation using LU for the case of two- 
qubit interactions. We recover the results of but 
through an alternative, simpler proof, based on known 
results of majorization theory — and thus avoiding the 
geometrical constructions of the original derivation 
The equivalence of LOCC and LU+anc strategies pre- 
sented in section HI, together with that of LU+anc and 
LU strategies for two-qubit Hamiltonians proved in [Q, 
imply that these results are also optimal in the context 
of LOCC, LO and LU+anc Hamiltonian simulation. 

We start by recalling some basic facts. Any two-qubit 
Hamiltonian H is equivalent, as far as LU simulation 
protocols are concerned, to its canonical form ^ 



where /ii > /12 > l^isl > and the operators cr^ are the 
Pauli matrices, 



(71 



1 

1 



, 0-2 



-i 

1 



,(73 



l-l)- (48) 



A brief justification for this canonical form is as follows. 
Any two-qubit Hamiltonian 



Ha® Ib + Ia® Hb + ^ 



(49) 



can efficiently simulate (or be efficiently simulated by) 
its canonical form (^): on the one hand we can always 
use traceless operators m and n as in Eq. ( ^0| ) to remove 
(or introduce) the local operators Ha and Hb; then the 
remaining operator hijai g) aj can be taken into the 
canonical form by means of one-qubit unitaries u and v 
such that (u (g) u) J2ij ^ij '^i ® ^3 {u^ ^v'^) is diagonal when 
expressed in terms of Pauli matrices. The coefficients hi 
in Eq. ( ^ ) turn out to be related to the singular values 
of the matrix hij . 

Therefore we only need to study the conditions for ef- 
ficient simulation between Hamiltonians which are in a 
canonical form. Let {I*!*!)} stand for the basis of maxi- 
mally entangled vectors of two qubits 



^-^(|oi) + |io)), |a>2) 



V2 



(|00>-|11», |<1>4) 



1 

V2 
1 

71 



(|00) + |11», 



(|01>-|10)).(50) 



Then H can be alternatively expressed as 



i7 = ^A,|$,)(ci>,|, 



(51) 



i=l 



where \i are decreasingly ordered, real coefficients fulfill- 
ing the constraint Ai = (coming from the fact that 
H has no trace) and 



Ai = /ii + /i2 - hz 

A2 = /li - /l2 + /l3 

A3 = -hi + /l2 + ft-s 

A4 = -hi - /l2 - ft.3- 



(52) 
(53) 
(54) 
(55) 



The most general simulation protocol using H and LU 
leads to 



H' = '^PkUk g) VkHul g) vl, 



(56) 



where we have assumed, without loss of generality, that 
H' is also in its canonical form, as in Eqs. ( ^7| ) and (js^), 
with corresponding coefficients ft-J and A^. 



Necessary and sufficient conditions for efficient 
simulation and optimal simulation factor 



^47^ Let us derive the necessary and sufficient conditions for 

H to be able to simulate H' using LU and for infinitesimal 
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simulation times. Uhlmann's theorem [ [LSl states that 
the eigenvalues of operator H' in Eq. (f36D , a unitary 
mixing of operator H , are majorized by the eigenvalues 
\i of H, that is 

A'l < Ai, 
A'l + A^ < A1+A2, 
A'l + A2 + A3 < Ai + A2 + A3, 
A'l + A2 + A^j + A^i = A1+A2 + A3 + A4, (57) 

where the last equation is trivially fulfilled due to the 
fact that H and H' are trace-less operators. Succinctly, 
we shall write A' -< A, as usual |19|. In terms of the 
coefficients hi and h'^ the previous conditions read 



K < hi, 
h[ + + /i3 < hi 



h2 - hs, 
h2 + hs, 



(58) 



and correspond to the s(pecial)-majorization relation, 
h' h, introduced in Ref. Q|. Thus, we have already 
recovered the necessary conditions of Q for H to be able 
to efficiently simulate H' in LU protocols pO| (and thus, 
since we are in the two-qubit case, also in LOCC proto- 
cols) . 

In order to see that conditions [and thus condi- 
tions ( |5^ ) ] are also sufficient for efficient LU simulation, 
we concatenate two other results of majorization theory. 
The first one (see theorem II. 1.10 of jl^) states that 
A' ^ A if and only if a doubly stochastic matrix m exists 
such that A^ = "Tiij^j- The second result is known 
as Birkhoff's theorem 19 , and states that the matrix m 
can always be written as a convex sum of permutation 
operators {Pk}, so that 



/a; 

Ai 
\ A^ 




(59) 



where each permutation affects two neighboring ele- 
ments. The corresponding three basic permutations of 
(<i>i, $2, ^4) can be easily obtained using LU. Indeed, 
in order to permute ($1, $2, ^3, ^4) into 



($2,3'l,$3,$4), 
($1,$3,$2,$4), 
($1,$2,*4,*3), 



(64) 



we can simply apply, respectively, the following local uni- 
taries: 



/ — iai 


/ 


— iai 






V2 


I + ia3 


/ 


- icr-i 




5 - 


V2 


I + i(7i 


/ 


— iui 


' 




V2 



(65) 



Therefore, any permutation a of the states (50) can be 
accomplished through local unitaries Tj, and any Hamil- 
tonian H' satisfying conditions ( ^8|) [equivalently, condi- 
tions (|5^)] can be efficiently simulated with H and LU. 

In the following we condense the previous findings into 
two results, Rl and R2, which provide an explicit answer 
to problems PI and P2, respectively, announced in sec- 
tion II. C of the paper. We assume that the two-qubit 
Hamiltonians H and H' are in their canonical form, with 
A, ft-. A' and h' the corresponding vectors of coefficients. 

Rl: Hamiltonian H' can be efficiently simulated by 
H and LOCC — or LU, LU+anc, or LO — if and only 
if conditions (j^j [or, equivalently, conditions (|5^/ are 
fulfilled, i.e. 

H' >Locc H ^ h' <,h ^ A' -< A. (66) 

R2: The simulation factor sh'\h for LOCC — or LU, 
LU+anc, or LO — protocols is given by the maximal s > 
such that sh' h or, equivalently, such that sX' -< A. 



This means that whenever conditions (|5^) are fulfilled 
we can obtain H' from H by using a mixing of unitary 
operations Ti, where each Ti permutes the vectors {!<&»)}, 



H' = Y,p,T,HT}. 



(60) 



Then, all we still need to see is that all 4! = 24 pos- 
sible permutations of the vectors {|$i)} can be per- 
formed through local unitaries Ti. Recall, however, 
that any permutation cr, taking elements (1,2,3,4) into 
(cr(l), (7(2), (t(3), cr(4)), can be obtained by composing 
(several times) the following three transpositions. 



(1,2,3,4)^(2,1,3,4), 
(1,2,3,4)^(1,3,2,4), 
(1,2,3,4)^(1,2,4,3), 



(61) 
(62) 
(63) 



B. Explicit optimal LU protocols 

The last question we address is how to actually con- 
struct a simulation protocol. That is, given H and H' , 
we show how to simulate sH' using H and LU, for any 

S g [Q,Sh'\h]- 

A complete answer to this question is given by a prob- 
ability distribution {pk} and a set of unitaries {uk ® Vk\ 
such that 



sH' 



E 

k 



PkUk ® VkHu\ ® v\, 



(67) 



where s S [Q,sh'\h\i and shi\h can be obtained using 
result R2. 

We already argued that it is always possible to choose 
all Uk ® Vk such that they permute the vectors of Eq. 
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(pO|), so that each Uk^Vk = is just a composition of 
the local unitaries of Eqs. (|65|). As before, let {Pfcj^li 
denote the 24 permutations implemented by the local 
unitaries {Tk}'^^i- Then the above problem reduces to 
finding an explicit probability distribution {pk\ such that 



k 

or, equivalently, such that 



(68) 



(69) 



This is done on appendix B using standard techniques 
of convex set theory. There we show how to construct a 
solution involving at most 4 terms pkTk for s < sh'\H: 
and at most 3 terms for optimal simulation, that is, when 
s = sh'\h- 

VI. CONCLUSIONS 

In this paper we have studied Hamiltonian simula- 
tion under the broader scope of LOCC protocols. We 
have focused on infinitesimal-time simulations, for which 
we have shown that LOCC protocols are equivalent to 
LU+anc protocols, also that LU+anc protocols are in 
general inequivalent to LU protocols (two-qubit Hamil- 
tonians being an exception). For two-qubit Hamiltoni- 
ans we have rederived and extended the results of Q , to 
finally provide the optimal solution using LOCC. 

Thus, the problem of simulating Hamiltonian evolu- 
tions has received a complete answer for infinitesimal 
times and using LOCC, for the simplest case of two-qubit 
systems. Several interesting questions remain open. On 
the one hand, the generalization of these results to sys- 
tems other than two qubits appears as challenging. On 
the other hand, the asymptotic scenario for Hamiltonian 
simulation, where H is used to simulate H' many times 
on different systems, certainly deserves a lot of attention. 

Finally, we note that entangled ancillary systems have 
been recently shown to be of interest in non-local Hamil- 
tonian simulation [|l). In particular, entanglement can 
act a catalyst for simulations, both in the infinitesimal- 
time and finite-time regimes, in that in the presence of 
entanglement better time efficiencies can be obtained, al- 
though the entanglement is not used up during the sim- 
ulation but is fully recovered after the manipulations. 
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APPENDIX A: EXTREME POINTS OF THE SET 
NON-LOCAL HAMILTONIAN SIMULATIONS 
ACHIEVABLE BY LU+ANC 

In this appendix we show that in LU+anc simulations 
any Hamiltonian of the form 

H' = {Oa'Ob'\U (E)V{H (g) Ia'b'W'' (E)V^\0a'0b') (A1) 

can be efficiently simulated by H, for any couple of uni- 
taries U and V acting on A A' and BB'. The result is 
valid also for more than two systems after a straightfor- 
ward generalization of the following proof. 

Notice that we can always write U and V using product 
basis {liAjA')} and {lisjs')} as 

= E E \'AJA'){^^J\ (A2) 

^ = E E Nsjs'X^yi, (A3) 

i=0 j=0 

where {If/iij)} and are other orthonormal basis of 

systems A A' and BB' , respectively, and denotes the 
dimension of system k. 

To perform this simulation, we need to make the out- 
put of the ancilla be the state \0a'0b'), unentangled with 
the systems AB. This can not be achieved by performing 
just transformations U and V, but by considering also a 
series of local unitaries {Ua ® Vb}, a G {0, • • • , dA' — 1}, 
b&{0,--- ,dB' - 1}, 



Ua = I®{J2 e^^'^lU'WA'DU, (A4) 

1=0 

Vt = I<E>{J2 e^^'-'llB'KlB'DV, (A5) 



and a constant probability distribution {pat}, Pat — 
l/{dA'dB')- Then we have that UI\Oa') = U^Oa'), and 
that J2a^a — dA'\OA'}{OA'\U , and similarly for Vb, so 
that we obtain 

Y^PabUa VbiH ® Ia'B'M ® V,^\0a'0b') = 
ab 

\0a'0b'){0a'0b'\U ® V{H ® Ia'B'W^ ® V^\0a'0b'){A6) 

Therefore Eq. ( [A^ ) defines a protocol that simulates the 
Hamiltonian of Eq. (Al) with unit time efficiency. 
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APPENDIX B: EXPLICIT TWO-QUBIT LU 
SIMULATION PROTOCOLS 



In this appendix we show how to find a probability 
distribution {pk} and permutations {Pk} such that 



A* = ^PkPk\ 



(Bl) 



for any two given four-dimensional, real vectors A and 
pi [pi ~ sX' in section V.B) such that fl ^ X, where 

We first note two facts that will allow us to use stan- 
dard techniques of convex set theory: (z) the set S = 
{f \ T ^ X} is convex, and (m) {PfeAjfl]^ are the extreme 
points of S, as it follows from Birkhoff's theorem [^ . 
We can then proceed as follows. 

Step (a): we check whether ft = PiX for any i — 
1, • • • , 24. If we find one such permutation we are done. 
Otherwise we move to step (6). 

Step (6): Facts (i) and (m) guarantee that there is at 
least one permutation Pfc, that we call Qi, and a positive 
e > such that 



fi = eQiX + (1 - e)f, 



(B2) 



where r also belongs to S, and therefore satisfies r -< 
A. In other words, we have to search until we find a 
permutation Qi such that 



iA)/(l-e)^ A, 



(B3) 



for some e > 0. Once we have found it we only need 
t o in crease e to its maximal value compatible with Eq. 
(B3). Let qi be this maximal value of e. Then we can 
write 



A* = 9i(3iA + (1 - qi)ii2, 



(B4) 



where p,2 < X \s on one of the surfaces of S — otherwise 
we could have taken a greater qi . 



Such a surface is, again, a (lower dimensional) convex 
set, whose extreme points are some of the PfeA's, and 
whose elements r fulfill f < X but with one of the ma- 
jorization inequalities replaced with an equality. This 
allows us to repeat points (a) and (6), but now aiming 
decomposing /i2 as a convex sum of vectors PfcA. That 
is, first we check whether fl2 corresponds to P^A for some 
k. And, if not, we search until we find a permutation P^, 
let us call it Q2, such that, again, 



(M2 - eQ2A)/(l -()<X. 



(B5) 



The maximum value of e compatible with this equation, 
say q, leads to a second term q2Q2 {12 = (1 — 91)9) for 
the decomposition of /i, 

M = + 92Q2)A-f (1 - gi - 92)^3, (B6) 

and to a new /i3, that lies on a surface of yet lower di- 
mensionality of the original convex set S. We iterate the 
procedure until the remaining vector pi lies on a convex 
surface of S of dimension zero, which means that the sur- 
face contains only one element, /2;. In this way we obtain 
the desired decomposition. 



y^gfcQfcA. 



(B7) 



k=l 



What is the minimal value of 11 For non-optimal simu- 
lation protocols we have that p — sA', where s < sh'\h, 
and p is in the interior of S, which is a three-dimensional 
set. Therefore the above procedure has to be iterated 
at most three times before we are left with a zero- 
dimensional surface of 5, and the minimal decomposition 
contains at most I = 4 terms. For optimal simulation 
protocols p = S[j,\fjX' is already in a surface of S, and 
therefore the minimal decomposition contains from 1 to 
3 terms. 
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